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We report on the fabrication and measurement of a graphene tunnel junction using hexagonal- 
boron nitride as a tunnel barrier between graphene and a metal gate. The tunneling behavior into 
graphene is altered by the interactions with phonons and the presence of disorder. We extract prop- 
erties of graphene and observe multiple phonon-enhanced tunneling thresholds. Finally, differences 
in the measured properties of two devices are used to shed light on mutually-contrasting previous 
results of scanning tunneling microscopy in graphene. 



The probability for an electron to tunnel through a 
very thin potential barrier depends strongly on the num- 
ber of available states in the target material at the same 
energy. Tunneling spectroscopy uses this principle to 
directly probe the density of states of complex mate- 
rials whose electronic properties are still poorly under- 
stood [1]. Furthermore, the phonon energy spectrum of 
a material can be determined by measuring the second 
derivative of the tunneling current [1]. Thus, a tunneling 
measurement can obtain information on both the elec- 
tronic and phononic properties of a material. 

In particular, tunneling experiments have been per- 
formed on graphene, a carbon monolayer with a honey- 
comb lattice where electrons are confined in two dimen- 
sions. In this material, electrons behave like relativistic 
particles obeying the Dirac equation, leading to a variety 
of exciting electronic behaviors [2] . In addition to prob- 
ing the density of states, prior tunneling spectroscopy 
measurements on graphene revealed features associated 
with interactions between electrons, phonons and disor- 
der [3-5, 7]. At low carrier densities, charge inhomo- 
geneities modify electron tunneling in scanning tunnel- 
ing microscopy (STM) [4, 5]. However, STM experiments 
have exhibited conflicting results. While some STM mea- 
surements have shown a strong suppression of tunneling 
at low energies [3, 6], attributed to interactions between 
electrons and K-point phonons in graphene [8], other 
experiments have not observed this effect, whether for 
graphene on Si02 [5, 9] or graphene on hexagonal-boron 
nitride (h-BN) [10]. 

In this Letter, we complement the existing tunneling 
experiments by fabricating graphene tunnel junctions us- 
ing h-BN as a tunneling barrier. We focus on two devices, 
referred to below as A and B. The devices were fabricated 
by the same process, but as we will see below they exhibit 
different tunnel behavior: one dominated by the density 
of states in graphene and one dominated by the energy 
dependence of the transmission probability. Each resem- 
bles a different set of prior STM experiments. Proper- 
ties of graphene extracted from the tunnel measurements, 
such as disorder-induced charge puddle size and Fermi 
velocity, agree well with previously reported results ob- 



tained via STM. In addition, we observe a rich electron- 
phonon interaction structure in the inelastic tunneling 
current, allowing for the identification of resonances near 
the phonon energies at the K and M-point in graphene 
as well as two low energy resonances that we attribute 
to Van Hove singularities in the h-BN phonon density 
of states. Finally, the energy dependence of the tunnel 
transmission is strongly influenced by the microscopic de- 
tails of the device geometry, shedding light on previously 
mutually-contrasting results in STM experiments. 

The use of h-BN as an insulating substrate has al- 
lowed for high-mobility graphene devices to be fabri- 
cated [11]. We exfoliate h-BN and graphene on sepa- 
rate degenerately-doped silicon wafer pieces capped with 
thermally-grown, 300nm-thick Si02- The thicknesses of 
the h-BN flakes used are measured with AFM to be 
s=s 2nm. After annealing both flakes in Ar and H2 at 
350 °C for 4 hours each with flow rates of 500 seem, 
the h-BN flakes are coated with a thick PMMA layer 
and the underlying SiC>2 is etched away with a solution 
of potassium hydroxide, which lifts off the PMMA with 
the h-BN attached to it. The PMMA-h-BN membrane 
is then attached to a micro-manipulator arm, similar to 
Ref. [11], and transferred on top of the graphene sheet. 
The PMMA is then dissolved away in acetone. Standard 
e-beam lithography techniques are employed to electri- 
cally contact the graphene with 10nm/50nm of Ti/Au 
and to add a top gate, which serves as a tunnel contact, 
on top of the h-BN flake [see inset of Fig. 1(a) for a 
schematic of the completed device] . 

The differential tunnel conductance g oc dl/dVr, 
where I is the tunnel current, is obtained as a function of 
the top-gate voltage Vr and backgate voltage Vq using a 
lock-in measurement with an excitation voltage of 1 mV 
at 13 Hz. Vq, applied to the degenerately-doped Si wafer, 
controls the average density of electrons n s in the sheet 
of graphene. Vt both modulates n s directly beneath the 
top gate and allows for the DC voltage bias dependence of 
the tunneling process to be probed. All experiments are 
performed at a temperature of 4.2 K. In-plane transport 
measurements on the devices show a charge-neutrality 
point (CNP) of 18(-22) V for device A(B). 
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FIG. 1: (a) (inset) schematic of the device structure used. 
An atomically-thin sheet of h-BN separates the top gate from 
graphene. An adjacent lead i is used to extract the current / that 
travels from the top gate (biased by Vt) to graphene. A voltage 
Vg is applied on the back gate to control the sheet density (n s ) in 
graphene. (main) The tunnel conductance g as a function of ■y/raj 
measured at Vt=0 and temperature of 4K. Dashed lines are guides 
to the eye. (b) g measured as a function of Vt for Vq=0. 



A plot of g as a function of y/n& for device A is shown in 
Fig. 1(a); here, n s is controlled by Vq, with Vr grounded. 
The tunnel conductance is proportional to \fn~ s away 
from the CNP, as shown by the dashed lines, directly 
reflecting the V-shaped behavior of the density of states. 
The variation of the density of states due to the metallic 
top gate is negligible within the range of Vr [1], and / is 
expressed as 



/(Vt) oc 



eVr 



p{E F + e)T(e, eV T )de 



(1) 



where p{E) is the density of states in graphene, Ep is the 
Fermi energy and T(e, eVr) is the energy and bias depen- 
dent electron tunneling transmission probability, which 
decays exponentially with distance and depends on the 
device geometry [12]. p(E) and Ep are given by: 



p(E) = 



2\E\ 



Tr(hvp) 2 
Ep = hvpkp = hvp^nris, 



(2a) 
(2b) 



where vp is the Fermi velocity of electrons. At Vr—0, the 
tunnel transmission is constant and g depends only on the 
density of states: g(Vr = 0, Vg) oc ep(Ep(VQ)) cx y/n^, 
as observed in Fig. 1(a). However, the density of states 
does not reach zero when n s vanishes but instead flattens 
out when n s ~ 5x 10 n c77i -2 , due to charged-impurity 
disorder on the graphene flake [13]. When n s becomes 
as low as the density of impurities, p or n-type charge 
puddles form around the impurities [14]. The slopes of g 
on the p-type and n-type regions differ by 20%, indicating 
that vp is 9% lower for electrons than for holes. This 
asymmetry is consistent with the presence of negatively- 
charged adsorbates on the sheet of graphene [15, 16]. 

The Vr-bias dependence of g for Vq=0 is shown in Fig. 
1(b). Results for device A show a smooth evolution of 
g with Vr, indicating no or relatively small variation in 
T(Vt), similar to Refs. [5, 9, 10]. In particular, there 
is no exponential suppression of g near zero-bias, as was 
observed in Refs. [3, 6, 7]. 

g(Vr, Vg) is moderately suppressed around Vr=0 and 
along a diagonal region [Fig. 2(a), outlined by white 
dashed lines]. This region corresponds to gate volt- 
ages such that n s is close to the CNP [18, 19], where 
the density of states is minimal. The tunnel probe 
also gates the graphene sheet and the relative back(top) 
gate capacitances Cq(Ct) can be determined by fitting 
CtVt + CqVq= constant. The slope of this line gives 
the capacitance ratio Ct/Cq ~ 72, which is lower than 
the ratio of 150 that one expects from a simple parallel 
plate capacitor model [12]. 

To observe the evolution of fine features in g, a deriva- 
tive with respect to the Vr-axis is taken, resulting in 
5g oc d 2 I/d 2 Vr [Fig. 2(b)]. We observe three collections 
of features in 5g for this range of voltages: peaks mov- 
ing along parabolas in Vt-Vg space away from the CNP 
(dashed curves), Vcj-independent peaks around Vr—0 
(dotted line) and strong features at negative values of 
Vt and Vg (squares), g for sample A is well approxi- 
mated by p(Ep — eVr) [12], therefore parabolic features 
in Sg(Vc, Vt) correspond to peaks in the density of states 
and follow contours of constant Ep — eVr • These fea- 
tures can be fit to determine vp (see [12] for informa- 
tion about the fit), obtaining 9.45 x 10 5 ± 8.5 x 10 4 m/s, 
in good agreement with the theoretically expected value 
of 1.1 x 10 6 m/s [17]. Two very strong diagonal peaks 
were obtained for Vr between ~-0.74V and -0.15 V in 
the range of Vg between -60 V and V (squares). These 
two peaks only appear in this device and are not sym- 
metric with respect to the CNP or Vr- While the origin 
of these two features is not known, it is likely that these 
are associated with a resonance in an impurity adjacent 
to the sheet of graphene. A simulation of the tunnel con- 
ductance in the presence of disorder is shown in inset of 
Fig. 2(b), reproducing the features indicated by dashed 
curves of Fig. 2(b) but not the features labeled by the 
dotted line or squares (see [12] for details of the simu- 
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FIG. 2: (a) <?(Vt, Vg) shows a suppression of g for the values of Vt 
and Vg outlined by white dashed lines, (b) (inset) Simulation of 
5<?(Vt,Vg). (main) Sg (Vt,Vg) measured for device A. Parabolic 
curves (black dashed lines) correspond to constant-density-of-states 
contours. In addition, a collection of VG-independent peaks are 
observed (one is highlighted by the dotted line). Squares indicate 
two unexplained features in Sg that are likely due to tunneling 
enhanced by an impurity adjacent to the graphene sheet, (c) Sg 
averaged over all positive values of Vq show the VG _m< Iependent 
peaks apparent in (a) and (b). Arrows correspond to the voltages 
10, 32, 68, 90 mV, from left to right. 



lation). The Vc-independent peaks, highlighted by the 
dotted line, are most clearly seen in Fig. 2(c) where a cut 
of Sg(Vr) that has been averaged over all positive values 
of Vg is shown. Four peaks are identified by arrows at 
10, 32, 68, and 90 mV (accurate to within lmV), with 
the former two peaks appearing stronger than the latter 
two. 

Similar measurements were performed on a second de- 
vice, B, where g primarily depends on Vt and varies only 
slowly with Vq , a result of the tunnel transmission vary- 
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FIG. 3: (a) Cut of <7(Vt) a t Vq=-10V for device B. showing a 
suppression of g for values of Vt less than ~70mV. (b) g(Vr, Vg) 
shows the threshold of 70 mV is independent of Vq for the range 
-50 V to 20 V. (c) Inelastic tunneling spectroscopy, measured by 
6g and averaged over the entire Vq range, shows peaks at 10, 40, 
66, 84 mV (shown by arrows for positive Vt but also apparent for 
negative values of Vt). (d) Peaks in Sg are VG _m< Iependent for the 
range -50 V to 20 V. 



ing faster than the density of states as a function of Vt- 
In Fig. 3(a), g is suppressed for |V"t| < 70 mV and a Vq 
dependence characteristic of graphene's density of states 
is only observed outside this suppression. The values of 
Vt where this suppression occured are independent of Vq 
for the range explored here [Fig. 3(b)]. A plot of Sg for 
device B [Fig. 3(c)] reveals peaks at similar voltages to 
those of device A, indicated by arrows at values of 10, 
40, 68 and 84 mV (±lmV), but with different intensi- 
ties than device A. These peaks, like for device A, are 
independent of Vg [Fig. 3(d)]. 

Features that are independent of n s and appear as 
peaks in d 2 I/d 2 Vr can only reflect the variations of the 
tunnel transmission and not the density of states, and are 
ascribed to phonon thresholds enhancing tunneling. One 
such threshold, at 68 mV for device A and B, has been 
previously measured and assigned to excitations of K- 
point phonons in graphene [3, 6]. The three other peaks 
are previously unreported in graphene tunneling experi- 
ments. The lower inelastic thresholds of 10 and 40 mV do 
not correspond to any phonon energy in graphite [20] and 
are likely scattering of tunneling electrons from phonons 
in h-BN. The phonon spectrum in h-BN exhibits a very 
flat ZA branch at 40meV along the MK direction [21]. 
This should be associated with a Van-Hove singularity 
(VHS) in the phonon density of states at 40meV, explain- 
ing why an inelastic threshold is observed at this energy. 
The 10 mV threshold is attributed to the ZA phonon at 
the A point [21], where a band flattening is also observed. 
The highest voltage kink (at 84 mV and 90 mV respec- 
tively on device A and B) is close in energy to both the 
M-point phonon in graphene [20] and another VHS in the 
h-BN phonon density of states [21]. We also observe that 
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FIG. 4: g for smaller values of Vt for (a) device A and (b) device B. 
Diamond-like features are attributed to charge puddle formations 
in the graphene sheet for each device. Coulomb blockade features 
are more pronounced in device B. 

peaks associated with h-BN phonons are much stronger 
for device A, whereas in device B, the graphene phonon 
thresholds are the strongest. 

The behavior of g at low Vt provides information on 
the role of impurities in the tunnel transmission (Fig. 4) . 
Both devices exhibit Coulomb diamonds. Similar fea- 
tures have been found in STM experiments [4, 5] and are 
due to the formation of charge puddles around charged 
impurities. These puddles behave like leaky quantum 
dots and dominate transport in graphene at low den- 
sity [13]. For device A, these Coulomb diamonds have 
different sizes, indicating that several charge puddles are 
involved. The number and size of these diamonds in- 
crease closer to the CNP, which is due to weaker screening 
of charged impurities [4]. The periodicity of the conduc- 
tance oscillations as a function of Vq for device A is used 
to extract the average capacitance of the dots to the back 
gate, resulting in a typical dot size ~ 200 nm 2 (see [12] for 
details of the dot size extraction) , in good agreement with 
values reported in Ref. [4]. For device B [Fig. 4(b)], the 
Coulomb blockade features are much more pronounced, 
indicating that only a few charge puddles are involved. A 
charging energy of 6 meV and a capacitance to the back- 
gate of 8.2 x 10 -21 F are measured, which corresponds 
to a dot size of approximately 70 nm 2 [12], smaller than 
that of device A. The ratio Ct/Cq for device B is found 



to be ~100 [12], larger than what was found for device 
A, which indicates that the tunnel gate is closer to the 
graphene sheet in device B. Given the exponential depen- 
dence of the tunnel current on the barrier thickness, one 
would expect the average value of g to be much higher 
for device B since the capacitance ratio is larger. This is 
not the case. This difference cannot be accounted for by 
the top gate geometry either, as the area of the top-gates 
for the two devices are comparable. This suggests that 
the tunneling area is much smaller for device B, indicat- 
ing that tunneling occurs primarily at a very localized 
point or collection of points. This localized tunneling 
can result from an impurity trapped between the h-BN 
and the graphene, locally enhancing the tunneling rate. 
Another scenario resulting in localized tunneling comes 
from a small imperfection in the h-BN lattice, causing 
the tunnel barrier to be lower in one small area. 

The data from device A and B show that the relative 
strength of elastic and inelastic tunneling strongly de- 
pends on the geometry of the tunnel junction since only 
B exhibits a strong suppression of g at small values of Vr 
[Figs. 3(a, b)]. To our knowledge, the phonon-enhanced 
tunneling effect has only been observed in STM exper- 
iments where the tips were prepared to be atomically 
sharp (see the supplementary infomation of Ref. [3]). 
This suggests that inelastic scattering plays a more im- 
portant role when the tip wavcfunction is spatially local- 
ized, and therefore has a very broad momentum distri- 
bution. h-BN enables us to observe tunneling across a 
two-dimensional interface: the tunneling electrons have 
a well-defined parallel momentum and inelastic tunneling 
is suppressed. 
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DEVICE GEOMETRY 

The graphene flakes are annealed in Ar/H2 at 350 °C 
and boron nitride flakes are transferred on top of them 
prior to any other processing, which allows for the in- 
terface between the two flakes to be very clean (See Fig. 
SI for an optical image of a completed device). The top- 
gated part of the graphene flake is several square microns 
large. However, the tunnel conductance is an exponen- 
tial function of the barrier thickness, so the effective tun- 
neling area depends strongly on the cleanliness of the 
interface. In fact, and as speculated in the main article, 
impurities in between the boron nitride and graphene can 
alter the nature of tunneling, resulting in the differences 
observed in the tunnel conductance g of device A and B. 

The measured capacitance ratio Ct/Cq ~ 72(100) for 
device A(B) is smaller than the theoretical value of 150 
given by a parallel plates model for this geometry. This 
difference is not fully understood but can be due to an 
effective dielectric constant lower than expected for the 
h-BN layer, or to an imperfect screening of the electric 
field of the gates by the graphene sheet. 



CHARGE PUDDLE SIZE 

Near the charge neutrality point in graphene, the den- 
sity of carriers breaks up into a series of n-and p-type 
puddles [1, 2], which behave as quantum dots [3]. The 
typical size of the charge puddles in our devices can be 
extracted from the Coulomb diamonds observed in Fig. 
4 of the main article. If one defines the capacitance ratios 
o-i = Ci/C total where i refers to the top gate (i=T) and 
back gate (i=G) capacitances, the edges of the Coulomb 
diamonds have slopes given by: 



dVo 1 — cct ^ —1 — «t 
dVr 2«g 2«g 



(1) 



From fits to the edges of the diamonds, it is therefore 
possible to extract the capacitance ratios Ct/Cq of 72 
and 100 respectively for device A and B. 

The capacitance of the puddle-induced quantum dots 
to the back gate is given by the periodicity of conduc- 
tance oscillations at Vr = and as a function of the 
back gate voltage. We substract a smooth background 



from g(Vr = 0, Vq) and calculate the fast Fourier trans- 
form of Sg. For example in the case of device A, we find 
a periodicity AVg of 7 V, which corresponds to a capaci- 
tance to the back gate of approximately 2 x 10~ 20 F. The 
back gate capacitance per unit area of the silicon oxide 
layer has been measured to be « 12 nF/cm~ 2 on different 
devices, which allows for an estimation of the typical dot 
size: 200 nm 2 in the case of device A, similar to that ob- 
served in Ref . [2] . A similar procedure is used to extract 
a puddle area of 72 nm 2 for device B. 



EXTRACTION OF THE FERMI VELOCITY 



The tunnel current can be expressed as: 



I(V T )oc [ p(E F + e)T(e,eV T )de. (2) 

J-eV T 



It follows that g is given by: 



gt = -^T « e P( E F - eV T )T{-eV T , eV T ) (3) 

r° d 

+ —p(E F + e)T(e,eV T )de 

J-eV T dV T 

Using the WKB approximation it is possible to esti- 
mate the tunnel transmission as a function of the barrier 
thickness d, the transverse electron mass in boron nitride 
to, the average barrier height U and the parallel momen- 
tum of the tunneling electron k//\ 



T(-eV T ,eV T ) = exp 



IdJlm 



U- 



(ft-Ay/) 2 

2to 




In our case, the barrier height U is comparable to half 
of the band-gap in boron nitride which we approximate 
by 4eV. Moreover, for electrons tunneling elastically at 
the K point, the parallel momentum k// is approximately 
equal to K w 1.7 A which is fairly high. As a conse- 
quence the tunneling energy eVt is small compared to 
the effective barrier height: U + ^'k//) and the tun- 
nel transmission T(—eVr,eVr) varies slowly with Vr as 
long as inelastic tunneling is negligible, as observed on 
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device A [4]. In that case, as a first approximation, one 
can neglect the variations of the tunnel transmission and 
write: 



9 = W t a pfyEp ~ eVr ^ 



dE F 
d(eV T ) 



-p(E F 



e)de 



1 d_ 
eV T de' 

P (E F - eV T ))(5) 



In graphene the Fermi energy is proportional to ^Jn 
and the derivative g^fc should diverge close to the 
charge neutrality point. However, we've seen that the 
density of states saturates at a constant value at low car- 
rier density [see Fig. 1(a) of the main article], and this 
divergence does not occur. As a consequence, the second 
term of Eq. 5 remains very small compared to p{E F — eV) 
for all applied voltages used in the experiment and the 
contours of constant tunnel conductance are very well 
approximated by curves of constant E F — eV. 

When the carrier density in the graphene sheet is large 
compared to the intrinsic doping no, the Fermi energy is 
given by: 



E F = hvpy/nn = hvp^J —(CtVt + CgVg + en ) (6) 

Using this expression it is possible to find an analytical 
expression for the contours of constant E F — eV. We find 
that these are parabolas of constant curvature: 



d 2 V G 



± 



2e A 



(7) 



Knowing the back gate capacitance, we can estimate 
the Fermi velocity from fits to these parabolas, and find 
a value of 9.45 x 10 5 m/s. 



The density of states is approximated by the theoreti- 
cal expression: pc{E) oc \E\ which is smoothly truncated 
to pg{E) = po under a cutoff energy that we take equal 
to 0.1 eV. Randomly placed Lorentzian peaks of random 
widths and heights are added to this expression to sim- 
ulate the resonant peaks we observed in Fig. 1(a) (main 
article). An example of the density of states we use is 
displayed on Fig. S2. We then numerically integrate this 
density of states to get the Fermi energy as a function of 
the carrier density, and calculate the tunnel conductance 
gt as a function of both gate voltages (Fig. S3) . We see 
that resonant peaks in the density of states give rise to 
two sets of curves: diagonal straight lines corresponding 
to constant Ep lines, and curves of constant Ep — eV, 
similar to the observed features of Fig. 2(b) (main arti- 
cle). 
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SIMULATION OF THE TUNNEL 
CONDUCTANCE 

Our goal in the simulation was not to determine the 
density of states for disordered graphene from ab-initio 
calculations, but to see what are the contours of con- 
stant tunnel conductance for a density of states resem- 
bling what was observed in Fig. 2 of the main article. 
To this end, we estimated the tunnel conductance from 
the WKB formula with an empirical density of states re- 
producing the main features in Fig. 2 (main article). We 
neglected variations in the tunneling transmission, as- 
suming that g oc pa(E F — eV t ). The Fermi energy itself 
is calculated by integration of the density of states. 




fig. Si: Graphene tunnel junction with a 2 nm thick h-BN tunnel barrier. The contacts and top gate 
are made by standard e-beam lithography (Ti/Au 10nm/50nm). The edges of the graphene flake 
are overlaid with dashed lines for clarity 




fig. S2: Simulation of the density of states in disordered graphene 
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fig. S3: Simulation of Gt as a function of the top gate voltage Vr and the back gate voltage V& s 



